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ABSTRACT 

We constrain the spectrum of the cosmic ultraviolet background radiation 
by fitting the observed abundance ratios carbon ions at 2 ~ 2-3 with those 
expected from different models of the background radiation. We use the recently 
calculated modulation of the background radiation between 3 and 4 Ryd due to 
resonant line absorption by intergalactic Hell, and determine the ratios of CHI 
to CIV expected at these redshifts, as functions of metallicity, gas density and 
temperature. Our analysis of the observed ratios shows that 'delayed reionization' 
models, which assume a large fraction of Hell at z ~ 3, is not favored by data. 
Our results suggest that Hell reionization was inhomogeneous, consistent with 
the predictions from recent simulations. 

Subject headings: cosmology: theory - diffuse radiation - intergalactic medium - 
quasars: general 
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Introduction 



Cosmological simulations of tlie intergalactic medium (IGM) have been used in the 
recent years to determine many important aspects of structure formation and reionization. 
One important input for these simulations, the cosmic ultraviolet background radiation, 
however, has remained uncertain. Our knowledge of the cosmic UV background, its 
intensity and spectrum and evolution with time, still lacks the precision that is needed to 
make the interpretations from cosmological simulations robust. 

It is believed that Hell reionization occurred much after the reionization of HI and 
Hel, with the help of hard UV photons from quasars, whose activity peaked around 
redshift ^ ^3. This is because of the relatively high io nization threshold of Hell (54.4 



eV) (e.g.. 



Madau fc Meiksin 



1994 : 



Sokasian et al. 



20021 ). Some recent observations may 



have 



(e.g. 



bund evidences for Hell Gunn-Peterson t rough in the spectr a of quasars at z > 2.8 



Jakobsen et al. 



of Hell reionization ( iMcOuinn et al 



1994; 



Schave et al. 



2000 



Smette et al. 



20021 ). Numerical simulations 



20091 ) also predict that the process is likely to be 
inhomogeneous, and accompanied by heating of the IGM. There are some evidences for 
such a heating from recent observations , although the interpretations are not yet c lear (e.g. 



Ricotti et al. 



2000 



Schaye et al. 



2000 



Theuns et al. 



2002 



Bernardi et al. 



20031). 



Recently, 



Madau fc Haardt 



(120091) pointed out an important aspect of the evolution of 
UV background radiation in the context of Hell reionization. They calculated the effect of 
resonant absorption by Hell Lyman series that is likely to significantly change the shape 
of the UV background radiation between 3 and 4 Ryd, whose magnitude depends on the 
abundance of Hell in the IGM. This modulation can be an important probe of the UV 
background radiation, when combined with the observations of absorption lines from metal 
ions such as CHI, whose ionization potential lies between energies corresponding to 3 and 4 
Ryd. 
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In this paper, we use the data from 



Agafonova et ah 



(120071 ) of column densities of CIII 



and CIV i n the redshift range of z ^ 2-3, to constrain the spectrum of the UV background 



radiation. 



Agafonova et al. 



(120071 ) used their observations to recover the shape of the 
spectrum with a Monte Carlo procedure, but their assumption for the shape of the spectrum 
did not include the possible attenuation between 3 and 4 Ryd. They claimed a sharp 
reduction in the flux of the UV background between 3 and 4 Ryd, which was interpreted as 
a sign of Hell Gunn-Peterson (CP) effect. In this regard, it is important to study the ionic 
ratios, especially those involving CIII, whose ionization potential falls between 3 and 4 Ryd, 
and the effect of bac kground radiation s of di fferent spectrum on these ratios. One of the 



predictions made by 



Madau fc Haardt 



(120091 ) was that the modulation of the spectrum by 



Hell absorption would significantly change the abundance ratio of CHI/CIV at 2; ~ 3, and 
a large CHI/CIV ratio could be explained without any need to invoke the Hell CP effect. 
We test these ideas in this paper. 



Calculations of ionic states 



We study the thermal and ionization evolution of gas exp osed to external ionizing 



radiation. Our calculation of the evolution of gas is similar to 



Gnat fc Sternberd (12007 



20091). A deta iled description of thermal and ionization evolution of gas can be found in 



VasilievI ( I2OIOI ). Here we briefly describe the method of calculation, present several tests of 



our code and discuss the choice of initial conditions. 



2.1. The description of the code 



A gas parcel is assumed to be optically thin for the ionizing radiation. We consider 
the time-dependent equations for all ionization states of H, He, C, N, O, Ne and Fe, 
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including all relevant atomic processes. Namely, we take into account the following major 
processes: photoionization, collisional ionization, radiative and dielectronic recombination 
as well as charge transfer in collisions with hydrogen and helium atoms and ion s. 
Atomic data for th e phot oionization cross section are adopted fr om 



Verner fc YakovlevI ( 119951 ). Auger effect probabilities taken from 



the recombination rate, 



Pequignot et al. 



fll99l[ ) 



br radiative recombination taken from 



Arnaud fc Ravmondl fll992[) and 



dielectronic recombination a dopted fr om 



ionization rate adopted from 



Mazzotta et al. 



Verner et al. 



Kaastra fc Mewe 



Verner fc Ferland ( 



Shull fc van Steenbergl (Il982[ ). for 



( 11998 



(11996 



19931) 



19961) 



), and for the collisional 



Voronovl (119971 ). for the ch arge transfer rates of ionizatio n and 



recombination with hydro gen and helium adopted from 



Kingdon fc Ferland 



1996|). 



Arnaud fc Rothenflug] ( 119851 ) and 



This system of time-dependent ionization state equations should be complemented by 
the temperature equation, which accounts for all relevant cooling and heating processes. 
Here we assume that a gas parcel cools isochorically, i.e. the cooling time is shorter than 
any dynamical time scale. Note that in the temperature equation we neglect the change in 
the number of particles in the system (which for a fully ionized plasma of hydrogen and 
helium remains approximately constant). 

The total cooling and heating rates ar e calculated using the photoionization code 



CLOUDY fver. 08.00, 



Ferland et al. 



19981 ). More specifically, we incorporate into the 



CLOUDY code a given set of all ionic fractions Xi calculated at temperature T, gas number 
density n and external ionization flux J(i^) and obtain the corresponding co oling and heatin 



rates. For the solar metallicity case we adopt the abundances re ported by 



( I2OO5I ). except Ne for which the enhanced abundance is adopted (IDrake &: Testa 



Asplund et al. 



20051 ). In 



all our calculations we assume the helium mass fraction to be Yue = 0.24. 



We solve a set of 68 coupled equations (67 for all ionization states of the elements H, 
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He, C, N, O, Ne and Fe and one for temperatur e) using a Variable-coefficient Ordinary 



Differential Equation solver flBrown et al 



19891 ). The time step is chosen to be the 



minimum value between {O.ltion, O.Oltcooi)' where tion and tcooi are the ionization and cooling 
times correspondingly. 

We study the ionization and thermal evolution of a lagrangian element of cooling 



gas exposed to extragala ctic spectra calculated 



cosmological simulations ( iMadau fc Haardt 



Dy using the radiative transfer in the 



20091 ) with and without the absorption in the 



He II resonant lines. Figure [U shows the extragalactic spectra used here. A brief description 
of the spectra is presented in Section 3. In our simulations we take into account the 
radiation spectrum from 1 to 10^ eV. 

The evolution of the gas certainly depends on a choice of the initial condition of gas, 
namely, the initial ionic composition and temperature. We consider such a dependence 
below (section 2.3). But ffist, we present several tests of the method of calculation. In these 
tests we start our calculations from temperature T = 10^ K with the collisional equilibrium 
ion abundances. During the calculation the gas number density is assumed to remain 
unchanged due to any dynamical or chemical processes. We stop our calculations when 
photoionization heating rate differs from cooling rate by less than 5% (|r — A|/A < 0.05). 
The ionic composition is almost "frozen" at this moment. Note that with such a stopping 
criterion a gas still remains "time-dependent" at the end of the computation, as it has not 
yet reached thermal and ionization equilibrium. 



2.2. Tests of the code 



Figures |2]|3] present the main tests of our code used in the calculations at the 
temperature range interested in this paper. Figure [2] shows the carbon ionization fractions 
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I-IV in the collisi onal equilibrium o 



original results by 



Mazzotta et al. 



jtained from our code (crosses of various types) and 



mm as well as the standard CIE test of the CLOUDY 
code. One can see a good coincidence between the results of our code and those obtained 
in the previous works. 

We show the comparison between the cases of photoionization equilibrium and 
time-dependent photoionization as well as collisional time- dependent and equilibrium with 
regard to GUI fraction in Figure |3l The case of photoionization equilibrium is shown with 
thick dash-dot lines and that of the time- dependent photoionization is shown by thick dot 
lines. We use a density of n = 10~^cm~'^ (lower thick lines) and n = 10~^cm~'^ (upper thick 
lines). The upper panels show the results for 10~^ Zq and those for the solar metallicity are 
shown in the lower panels. In the photoionization calculations the extra galactic ionizing 



spect rum without any modulation - HM model, (see Figure [T]) at z = 3 (IMadau fc Haardt 
20091 ) is assumed. 



Figure [3] also shows the collisional tim e- dependent GUI fraction obtained using our 



code (thick solid line) and that taken from 



Gnat fc Sternbergl (120071 ) (thin sohd hne), as 



well as the GIE fractions obtaine d from GLOUDY ( t hin da sh line). Firstly, one notes the 



difference of our results from the 



Gnat fc Sternberg] (120071 ) data. This difference can be 



explained by the fact that we use a different set of atomic data (radiative and dielectronic 
recombination rates) for several ions. We use the data for radiative recombination from 
Verner & Ferland (1996), which partially included the refitted data from previous works 
(see a b ove for full list of papers), and for dielectronic recombination from 



references in 



Mazzotta et al. 



Gnat fc Sternberg 



(119981 ). iGnat &: Sternbergl (120071) used more recent data of Badnell and coauthors (see 



2007 



Bryans et al. 



20061 ). However, our GIE results are 



close to those obtained in GLOUDY code, which traditionally includes the most recent 
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atomic data^. So we assume that such difference may be mainly explained by small changes 
in thermal evolution. For example, it is well known that the dielectronic recombination 
depends exponentially on temperature, so we expect that for higher temperature and lower 
metallicity this difference becomes smaller (see Figure [3]). 

Secondly, we should note that both equilibrium and time-dependent collisional ionic 
composition has no dependence on gas density, because all processes are two-body. In 
contrast, the ionic composition in the photoionization case strongly depends on gas 
density. Moreover, the ionizing radiation forces the ionic composition of gas to settle 
on to equilibrium. In low density gas the ionic composition in the time- dependent 
photoionization case is expected to differ strongly from that in collisional both equilibrium 
and time-dependent cases, but it tends to the ionic fractions in the photoequilibrium. In 
high density gas we expect that the time-dependent photoionization ionic composition 
tends to the time-dependent collisional one, whereas the photoequilibrium ionic fractions 
should be close to those in the CIE. In addition, the ionic composition strongly depends 
on metallicity. For low metallicity the difference between time- dependent and equilibrium 
ionic fractions is expected to be small, but it increases with metallicity. 

In Figure [3] one can see that the GUI fraction in the time-dependent photoionization 
model does not fully coincide in the whole temperature range with either the photoionization 
equilibrium or pure collisional models. For example, at low metallicity, Z = 10^'^ Zq, for 
high density value, n = 10~^ cm~^, the time- dependent photoionization GUI fraction shows 
some difference from the photo equilibrium case. In contrast, for n = 10~^ cm~^, the GUI 
fractions in the time-dependent photoionization and photo-equilibrium cases are close to 
each other in the whole temperature range. For solar metallicity, the difference between 
time-dependent and equilibrium photoionization becomes more significant for both number 



http: / / wiki.nublado.org/wiki /RevisionHistory 



- 9 - 



density values. Also one should note from the right panels of Figure [3] there is almost no 
dependence on metallicity in equilibrium photoionization case. The difference between the 
CIII fraction for 10^'^ Zq and solar metallicities is only in the value of temperature reached 
in the cooling process: for higher cooling rate (solar metallicity) this temperature is lower. 

As it is expected for higher density value, n = 10~^cm~^, the CIII photoionization 
fractions are close to the coUisional ones in both equilibrium and time-dependent cases in 
a wide temperature range. For example, the ionic fractions in time-dependent coUisional 
(thick solid line) and time- dependent photoionization cases (thick dotted line) almost 
coincide for logT ^ 4.7 at Z = 10~^Zq and for logT ^ 4.18 at solar metallicity. In 
summary, both time-dependent coUisional and photoequilibrium in the presence of a 
significant ionizing radiation flux like the extragalactic background models produce the 
ionic composition, which differs from that in time-dependent photoionization model. Due to 
the very non-linear dependence of ionic fractions on temperature, metallicity and density it 
is difficult to say where it is possible to use time-dependent coUisional or photoequilibrium 
models, so in the present work we use the more complex, but more adequate time-dependent 
photoionization model. 



2.3. Dependence on the initial temperature 



Here we should consider a choice of the initial temperature and ionic composition for 

" t range z = 2 — 3. M any cosmological simulations 



our study of 


absor 


□ers at the redshi 


ft 


(Dave et al. 


2001 




Cen & Fane 


2006 





Bertone et al. 



20081 ) show that the intergalactic gas 
at z = 2 — 3 mainly resides in the diffuse phase with a typical temperature T ^ lO^K and 
overdensity 6 < 1000, and a small part of gas (in the best model this fraction is less than 
5%, see e.g. Dave et al. 2001) can be found in the warm- hot phase with lO^K ^ T ^ lO^K. 
It is expected that a gas enriched by metals and expelled from galaxies would have passed 
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through strong shock waves. Since the typical velocity of galactic winds varies from several 
tens to hundreds km s~^, the intergalactic gas enriched by metals is likely to go through 
such shock waves. We assume that such a gas was initially heated by shocks up to the 
temperature T ^ lO'^K. We are interested in subsequent evolution of this shock-heated gas. 
Here also we should note, firstly, that the heating rate produced by the strong ionizing 
background at 2; = 2 — 3 does not allow low-density gas with n ~ 10~^cm~^ to cool 
effectively below T ~ 3.5 x lO^K, and secondly, that here we do not include adiabatic 
cooling due to the expansion of the universe, which dominates over the radiative cooling for 
very low density gas, n ~ 10~^ cm"'^ at 2; = 3. 

In the following set of calculations we study the dependence of the ionic composition 
evolution of a gas on the initial temperature Tj. We vary the initial temperature in a wide 
range Tj = 5 x 10"^ — 10®K. Here we follow the evolution of gas irrespective of the time 
needed to cool from Ti to the temperature value, when our stopping criterion (mentioned 
above) is reached. Further we compare the cooling time with the comoving Hubble time 
at 2; = 2 — 3. The gas cooling is strongly determined by both initial conditions of the gas 
(density, metallicity and temperature) and the UV background radiation. We assume the 
initial ionic composition corresponds to Tq = 2 x lO^K. Taking such conditions we simulate 
a gas parcel with Tq = 2 x lO^K that initially passes through a strong shock wave front 
with Tg ^ Tq. Because of the short ionization timescale of metals (ionization fraction of 
hydrogen reaches almost unity for Tq = 2 x lO^K in both CIE and time-dependent cases) 
the ionic composition of a gas passed through a shock front with temperature ^ Tq is 
expected to tend to that with Tj = lO^K at T ^ Tg. Certainly, if the initial temperature is 
low enough, then the ionic composition should strongly depend on the initial temperature 
value. 

Figure m presents the CHI fraction [upper panels) and the CHI/CIV ratio (lower 
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panels) for different starting temperature values and the ionization composition initially 
corresponding to that in the CIE at Tq = 2 x lO^K. In the upper panels (the left is for 
Z = 10~^Zq, the right is for the solar metallicity) the significant deviations from the GUI 
fraction for the evolutionary track with Tj = lO^K is found for T < lO^K, whereas the GUI 
fraction for tracks with Tj = lO^K and Tj = lO^K almost coincides. The time needed for 
cooling a gas from Tj ~ lO^K with n ~ 10~^ cm~^ and Z ~ IO^^Zq is comparable with 
the comoving Hubble time at 2; = 2 [tniz = 2) ~ 1.1 x 10^^ s). The cooling times from 
Ti ~ lO^K is slightly shorter than the Hubble time at z = 3, but the GUI fraction strongly 
differs from that for Tj = lO^K. But here we are interested in the ionic ratios, and not in 
column densities of individual ionic species. 

In the lower panels of Figure H] the GIII/GIV ratio demonstrates a weak dependence on 
the initial temperature. The ratio for Tj = 4 x lO^K almost coincides with that for higher 
temperature values. For the lowest temperature value, Tj = 5 x lO^K, considered in this set 
of calculations, the ionic ratio settles to the common trend almost at the same temperature 
T 5 X lO^K for n = 10~^ cm~^ (in this case at the beginning the photoheating is 
significant) and at T ^ 3.2 x lO^K for n = 10~^ cm~^. One can see the vertical parts of the 
GIII/GIV ratio tracks. However, the time needed for settling on the common trend (the 
time at the vertical part of the track) is significantly lower than the time needed for reaching 
the stopping criterion. For instance, for a gas with n = 10"^ cm~^ and Z = 10~^Zq such a 
timescale is about 2 x lO^^s in case of the initial temperature = 5 x lO^K, whereas the 
time of the calculation before the stopping criterion reached is more than order greater, 
~ 5.2 X lO^^s^ O.Ttfjiz = 3). The timescales for higher density or metalli city are smaller 



Thus , for T ^ 4 X lO^K (that is consistent with the line widths inferred by 



Agafonova et al. 



20071 ) the GIII/GIV ratio in a gas with initial temperature Tj ^ 5 x lO^K almost coincides 
with that for Tj = 10^ — lO^K, and the timescale of calculation needed for reaching our 
stopping criterion is smaller than the comoving Hubble time at z = 2 — 3. In light of this, we 
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can start our calculations from temperature T = 10^ K and consider the obtained CIII/CIV 
ratio. We should emphasize that we are interested in the temperature range T ^ 4 x lO^K, 
where the ionic ratios for Tj = 5 x 10^ — lO^K are very close (see lower panels of Figure H]), 
thus, our conclusions are almost independent on the initial temperature value, if this value 
is greater than 5 x lO^K. Further analysis of the dependence of ionic composition on the 
initial conditions is out of scope of this paper and will be done elsewhere. 



3. CIII/CIV ratio 



In this section we study the influence of the extragalactic ionizing background with and 
without Hell absorption on Clll/ClV ratio. Nonequilibrium ionization states are calculated 
for ultraviolet (UV) spectra at redshifts z = 1.87,2.48,2.65,2. 83,3 (all spectral data w ere 



Madau fc Haardt 



fl2009[ ). 



kindly provided by F. Haardt) for three models presented in 
hereafter MH09 (the abbreviations of models are the same as in MH09): a) absorption in 
the Hell resonant lines was neglected - HM model, b) the sawtooth modulation was added 
- HM + S model, where HeII/HI=35 in optically thin absorbers, and c) the model where 
the ratio Hell/Hl were artificially increased to 250 - DR (delayed reionization) model (for 
details see MH09). 



Note that the MH09 spectra differ from the iHaardt fc Madaul (119961 . |200l|) spectra. 
First, MH09 spectra include quasar contribu tion only, and, second , in the UV background 



calculations a QSO luminosity function from 



Hopkins et al. 



( 120071 ) is used, which produces 



a very steep decline of the ionizing flux at high redshifts. Figure [T] shows the far UV part of 
the MH09 spectra for z = 2.48. 

The left panel of Figure [5] shows Clll/ClV ratios in the equilibrium calculation (using 
CLOUDY) for Z = 10~^Zq (but there is no dependence on metallicity in the equilibrium. 
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see Section 2.2 and right panels of Figure |3]). The two middle panels of Figure [5] shows 
CIII/CIV ratios in the nonequilibrium calculation (using our program for nonequilibrium 
calculation) for Z = 10~^Zq and solar metallicities, correspodingly, for UV background 
spectrum aX z = 2.48. The difference between models without (HM) and with {HM + S) 
the inclusion of the sawtooth modulation is small for the whole range of metallicities and 
densities. The ratios for the HM and HM + S models differ from each other only by a 
factor of ~ 1.3 for low density at T ^ 10^ K and they almost concide for n = 10~^cm~^ 
or at T > 10^ K. In the DR model the ratio A^(CIII) /iV(CIV) is greater than those in the 
HM and HM + S models by a factor of ~ 10 at T ^ 10^ K. This factor decreases with the 
increase of density and metallicity. For example, for n = 10~^cm~^ a significant difference 
can be found only at T ^ 6 x 10^ K for Z = 10~^Zq and T ^ 2 x 10^ K for solar metallicity. 

For high density {n ^ 10~^ cm~'^), the ratio is close to that in the coUisional limit at 
T ^ 5 X 10^ K. One can therefore conclude that large sawtooth modulation as present in 
the DR model leads to higher abundances of GUI ions by a factor of 10-1.3 for gas density 
n = 10~^ — 10~^ cm~^, respectively. The increase of GUI and GIV abundances in the DR 
model arises from the decrease of GV and GVI abundances, because of lower ionization flux 
in the energy range ~ 50 — 500 eV in comparison with two other models, so the growth of 
the GUI abundance is higher than that of GIV making GIII/GIV ratio higher. 

Our results firstly show that the ratio of GIII/GIV is not s i gnific antly changed by the 



sawtooth modulation of Hell, as supposed by 



Madau fc Haardt 



( 120091 ). Since the ionization 



threshold of GUI (47.9 eV) lies within the Ly/3 absorption feature from Hell, and that of 
GIV lies well beyond the range where Hell resonant absorption changes the spectrum, it 
was expected that observations of the GIII/GIV ratio along with the theoretically calculated 
spectrum would be a good prob e of the physical conditio ns in the IGM at z ~ 3. Our 



detailed calculations confirm the 



Madau fc Haardt 



(j2009[ ) claim that the ionization rate of 
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CIII does not significantly differ between the "HM" and "HM+S" spectra, although the 
"DR" spectrum does make a difference. In other words, the CIII/CIV ratio is a good probe 
only for distinguishing between the extreme cases of " HM" or " DR" , and not for probing 
the attenuation caused by standard abundances of Hell in the IGM. 



4. Statistical analysis 



Agafonova et al. I fl2007l ) give measurements of CIII and CIV column densities for 10 
absorbers in the redshift range 2 < z < 3. We use the ratio R = A^(CIII)/A^(CIV) to 
constrain the ionizing radiation at far-UV range (Table 1). 

For our analysis we consider three models of ionizing radiation, as discussed in the 
previous section, in the redshift range 2 < z < 3 , the metallicity in the range Z = 10~'^-1 
the number density in the range n = 10~^-10~^ cm~^, and temper ature T < 4 x 10^ K. Thi s 
range of temperature is consistent with the line widths inferred by 



Agafonova et al. 



((20071). 



The measured ratio -Rob depends on the parameters n, T, Z and the model of ionization. 
One expects the Z, n and T to vary appreciably from one absorber to another and therefore 
the usual approach cannot be applied. Only the background ionizing flux can be assumed 
to vary slowly with redshift. In the light of these assumptions, we define, for each absorber, 
e = l-Rth — Rob\/\^R\, here i?th are the CIII/CIV ratios from our models, i?ob are the the 
observed values, and AR are the corresponding errors on these measurements (Table 1). 
For each absorber, we search the entire parameter space of T, n, Z, for a given model of 
ionizing radiation corresponding to the redshift of the absorber. For any absorber, a model 
is deemed acceptable if e < y/E (that corresponds to Sa- level of the theoretical model). 

Our results for the three models are shown in Tables 2-4. As seen from the tables, t he 
models HM and HM + S can fit well for 8 out the 10 data points of Agafonova et al. ] 
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( 120071 ) ■ However, the DR model is acceptable for only six out of the ten absorbers. One 



absorber is fit well by none of the models (absorber 4). 



In Tables 2-4 we also show the results for the equilibrium models using CLOUDY 
(Figure 3). The results in the two cases have qualitative similarities, e.g. absorber 4 is ruled 
out by all models. However, there are also quantitative differences. Many of the absorbers 
ruled out by the nonequilibrium models are now allowed or vice versa. An interesting 
departure between the two results is that the DR model is ruled out for as many as 8 
absorbers. 

Next we briefly discuss how our results depend on the choice of n, T, and Z. The 
preferred range of temperature is generally > a few x 10'^ K to obtain acceptable solutions 
for all the absorbers. We do not list the best fit values of metallicity in the tables as the 
metallicities in the entire range we considered are allowed. We also note that this is also 
expected if the CHI/ CIV ratios were obtained from equilibrium ionization models. We 
note that the dependence on n is stronger; the value of n is greater than 10~'^cm~'^ for 
all the acceptable models, except for the DR model. This is in line with the fact that 
the absorbers have densities comparable to or larger than virialized haloes at 2 ~ 3, as 
distinct from smaller column density Lyman- o; clouds which are only mildly overdense. 



Agafonova et al. 



(120071 ) found a significant HI column densities. 



In the measurements 

iV(HI) ~ 10^^~^^cm~^. Also in the majority of the absorbers they detected singly- ionized 
carbon and silicon: A^(CII) ~ lO^^cm"^, A^(Sill) ~ lO^^cm"^, whose ionization potentials 
are below 13.6 eV. But detection of singly-ionized or neutral (like C, O or Si) species 
depends on the filling factor of cold clouds in the ISM of absorber and therefore it is not 
guaranteed that such a cold cloud must appear along the line of sight. 

Even though this evidence is not conclusive to rule out the DR model, it is probably an 
indication that this line of enquiry is likely to better constrain the background far-UV flux. 
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We note that in absorbing systems wi th iV(HI) > 10 



16-17 



radiation above 54 eV may be attenuated (IHaardt fc Madau 



cm ^, the UV background 



1996 



McQuinn et al. 



20091). 



That can lead to higher ionic ratio N(CIII)/N(CIV), e.g. for the absorber no. 10 (see 
Table 1), thus the absorber no. 10 may suffer from self-shielding. But it is unlikely that 
the UV flux is severely attenuated in the case of other absorbers, because in the extreme 
case of zero flux, the ionic ratio N(CIII)/N(CIV) is expected to be larger than ~ 20 for 
T < 4 X lO^K, according to Fi gure[5](the thin dot-dashed line for the coUisional limit), and 
all absorbers except no. 10 show a smaller ratio than this, suggesting that attenuation of 
UV flux is unlikely to be a major problem. 



5. Discussions 

Our results suggest that it is possible to explain the ionic ratio without resorting to 
additional attenuatio n of UV flux between 3 and 4 Ryd by continuum Hell absorption, 
or a Hell GP effect. 



Agafonova et al. 



(120071 ) used their data to reconstruct the far-UV 



background radiation spectrum, and claimed a value of rGp(Hell) ~ 2.5-3. This was 



compared with other direct obse rvations of Hel 



opacity rGp(Hell) ~ 2.09 ± 0.1 (IReimers et al. 



GP a bsorption trough, at z = 2.87, with 



20051 ). It is difficult to disentangle the 



effects of absorption from a continuous IGM and an ensemble of clouds, as in the case of HI 
GP effect (see 



Becker et al. 



20071 ). although 



McQuinn 



(120091 ) has argued that the case for 



Hell GP effect is stronger than this, owing to the low abundance of Hell compared to HI. In 
light of the difficulty of identifying true GP effect, our results of explaining the ionic ratios 
of CHI/CIV with standard attenuation from Hell absorption suggest that the observations 
(120071 ) do not necessarily need a Hell GP effect, although it cannot be 



of 



Agafonova et al. 



completely ruled out. 



Our results also clearly show that it is difficult to select one model of far-UV flux from 
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the observed data. It is possible that the measured ratios are sensitive to the local conditions 
inside the cloud. It is also possible that the far-UV radiation field is inhomogeneous over 



spheres correspondi ng to this species might jus t be merging at 2 : 3 



Smette et al. 



McQuinn et al. 



~ 3 


Jakobsen et al. 


1994; 


(2009 


) have also shown that 



2OO2I ) . Numerical simulations by 
the reionization of Hell remains patchy even at 2 ~ 3. Therefore, it is entirely conceivable 
that it is inhomogeneous reionization and not fiuctuation in the density, metallicity and 
temperature that is responsible for the wide range of observed CIII/CIV ratios. 

A bigger data set involving different ionic ratios is however needed to address some of 
these questions in detail, and confirm the patchiness of the Hell reionization process. Also, 
to put better constraints on the spectrum profile we need to study ionization structure of a 
cloud coupled with gas dynamics and radiation transfer. Radiation transfer effects would 
particularly affect large HI column density clouds, and change the predictions of ionic ratios 
in them, but these effects are outside the scope of the present paper. Furthermore, the 
predicted ionizing background depends on the spectrum and luminosity function of sources, 
and this introduces an additional uncertainty in analysis of the kind we have presented here. 
It is possible that the effect of these uncertainties would be mitigated by studying different 
ionic ratios, and we will report the advantages of using other ions in a future paper. 



6. Conclusions 



In this paper we used the data of 



Agafonova et al. 



(120071 ) for abundance ratio of CHI 



to CIV in quasar absorption systems at z ~ 2.4-2.9, to constrain the sp ectrum of the cosm ic 



UV b ackground radiation. We have used the spectrum calculated by 



Madau fc Haardt 



2OO9I ) which took into account the modulation from Hell Lyman series lines, and also 



models which assume a large abundance of Hell by artificially increasing the Hell/HI ratio. 
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Our results can be summarized as follows: 

1. The ratio between CIII/CIV is not a sensitive probe of the attenuation of the far-UV 
spectrum between 3 and 4 Ryd for standard abundance of Hell in the IGM. 

2. The observed data do not favor the models with additional abundance of Hell in the 
IGM (the "DR", delayed reionization models). 

3. There is a large variation in the fitted values of density and temperature for the 
absorption systems. This is indicative of an inhomogeneous reionization of HeH at 
these epochs, that has been suggested from numerical simulations, and our results 
lend support to these models. 
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Fig- !•- 



The 



(iMadau fc Haardt 



ar u ltraviolet part of the extragalactic ionizing spectra at z 



20091 ): the spectrum without absorption in the Hell resonant lines - 
HM model - thick dot-dashed line, the spectrum with the sawtooth modulation - HM + S 
model (HeII/HI=35 in optically thin absorbers) - thick dashed line, and the spectrum model 



where the ratio Hell/HI were artificial. 
- thick solid line. For comparison the 



y increased to 250 - Df ? (dela yed reionization) model 



Haardt fc MadaiJ (11996 



20011) spectrum {EM^V) for 



quasar-only (thin dot-dashed line) is added. 
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this work X 
Mazzotta et al 




4- 4.2 4.4 4.6 4.8 5 



log T, K 



Fig. 2. — The coUisional equilibrium carbon I-IV state frac tions for our code (crosses) are 
superposed on the original results by iMazzotta et al. I ( 1l998l ) (solid lines) and the standard 
collisional equilibrium test of CLOUDY (dash-dotted lines). 
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this work . 
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4.5 5 

log T, K 



4.5 6 
log T, K 



Fig. 3. — The CIII ion fraction for 10 ^Zq (upper panels) and solar (lower panels) metal- 
licities for coUisional (left panels) and photoionization (right panels) cases. The collisional 



time-depen dent ion fraction o 



taken from 



j tained using our code is shown by thick solid line and that 



Gnat fc Sternberg 



(120071 ) is depicted by thin solid line (GS07). In both left 
panels the CIE fraction of CIII obtained from CLOUDY is shown by thin dashed line. The 
time-dependent photoionization (thick dot lines, "photo-NEQ") and photo equilibrium (thick 
dash-dot lines, "photo-EQ") fractions are shown for n = 10~^ cm~^ by lower thick lines and 
for n = 10~^ cm~^ by upper thick lines, correspondingly. In the photoionization calculations 
the ex tragalactic ionizing spe ctrum without any modulation {HM model in Figure [1]) at 



z = 3 ( iMadau fc Haardt 



20091 ) is assumed. 



4.0 4.2 4.4 4.6 4.8 5.0 5.2 4.0 4.2 4.4 4.6 4.8 5.0 5.2 
logT, K logT, K 



Fig. 4. — Upper panels. The CIII ion fraction for gas with Z = IO^'^Zq (left panel) 
and with solar metallicity (right panel), and n = 10~^ cm~^ (solid lines), n = 10~^ cni"'^ 
(dashed lines) for different initial temperature values. Thick lines correspond to the initial 
temperature T = lO^K, thin lines present the evolutionary tracks for initial temperature 
T = 5 X 10^, 10^ 4 X 10^ lO^K (from top/left to bottom/right). Lower panels. The CIII/CIV 
ratio for gas with the same parameters as in the upper panels. 
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2.4 2.7 3.0 



Madau fc Haardt 



fl2009h 



Fig. 5. — Left panel. The equilibrium CIII/CIV ratio for the 
spectrum aX z = 2.48 for Z = IQ^^Zq metalhcity is plotted as a function of gas temperature 
for various cases. The solid lines depict the ratio for the "HM" extragalactic spectrum, 
dashed lines - "HM+S", and dotted lines - "DR" models for three values of number density 
n = 10~^, 10~^, 10~^ cm~^ (from bottom to top). For clarity, the models with n = 
10~^ cm~^ are shown by gray lines. Thin dot-dash line corresponds to the ratio in the 
coUisional limit: for Z = 10~^Zq. Left middle panel. The same as in the left panel, but for 
nonequilibrum models. Right middle panel. The same as in the left middle panel, but for 



solar metallicity. Right panel. The observational data of CHI and CIV co^ 



10 absorbers in the redshift range 2 < 2; < 3 as given by 



Agafonova et al. 



umn densities for 



fl2007h . 



Table 1. Carbon ionization ratios for absorption systems 



absorber Zabs R = Ncm/Ncw error in R 



1 


2.379 


1.438 


0.160 


2 


2.433 


6.944 


0.449 


3 


2.438 


2.680 


0.082 


4 


2.568 


1.000 


0.083 


5 


2.735 


3.083 


0.076 


6 


2.739 


15.681 


1.782 


7 


2.741 


3.158 


0.499 


8 


2.875 


4.750 


0.025 


9 


2.939 


16.363 


2.645 


10 


2.944 


21.765 


0.311 
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Table 2. The best parameters for the model HM. Columns 2-4 present results for the 
nonequilibrium model. The next three columns show the results for the equilibrium model 



absorber T (K) n (cm ^) e T (K) n (cm ^) e 



1 


40000 


0.001 


4 


31622 


0.001 


2.1 


2 


27200 


0.01 


0.1 


39810 


0.001 


8.2 


3 


30300 


0.001 


0.073 


31622 


0.001 


1.82 


4 


21740 


0.0001 


6.9 


31622 


0.001 


9.2 


5 


35460 


0.001 


0.02 


39810 


0.001 


2.5 


6 


33300 


0.01 


0.027 


19952 


0.01 


0.18 


7 


36370 


0.001 


0.005 


39810 


0.001 


0.24 


8 


37920 


0.01 


1.36 


39810 


0.001 


58 


9 


12840 


0.01 


0.006 


25118 


0.01 


0.2 


10 


31730 


0.01 


0.15 


25118 


0.01 


1.5 
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Table 3. Same as Table 2 for the model HM+S 



absorber T (K) n (cm e T (K) n (cm ^) e 



1 


39870 


0.001 


5.7 


31622 


0.001 


2.83 


2 


28210 


0.01 


0.12 


39810 


0.001 


6.68 


3 


34930 


0.001 


0.07 


25118 


0.001 


0.85 


4 


23800 


0.0001 


6.1 


31622 


0.001 


10.7 


5 


30200 


0.001 


0.1 


39810 


0.001 


2.29 


6 


36310 


0.01 


0.04 


19952 


0.01 


0.4 


7 


29410 


0.001 


0.02 


39810 


0.001 


0.37 


8 


23400 


0.001 


0.4 


39180 


0.001 


31 


9 


13670 


0.01 


0.032 


19552 


0.01 


0.012 


10 


37070 


0.01 


0.015 


39810 


0.01 


0.12 
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Table 4. Same at Table 2 for the model DR 



absorber T (K) n (cm ^) e T (K) n (cm e 



1 


39810 


0.0001 


0.24 


31622 


0.001 


42 


2 


31100 


0.001 


0.02 


31622 


0.001 


3.1 


3 


16840 


0.0001 


3.8 


31622 


0.001 


68 


4 


39810 


0.0001 


5.7 


31622 


0.001 


87 


5 


16840 


0.0001 


9.4 


31622 


0.001 


68 


6 


37830 


0.001 


0.01 


39810 


0.001 


0.14 


7 


16840 


0.0001 


1.6 


31622 


0.001 


10 


8 


39990 


0.001 


10.8 


31622 


0.001 


142 


9 


16710 


0.01 


0.001 


39810 


0.001 


0.16 


10 


39630 


0.01 


0.015 


39810 


0.001 
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